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Abstract. - In the first part we discuss how the BEC picture for magnons is modified by 
anisotropics induced by spin-orbit coupling. In particular we focus on the effects of antisym- 
metric spin interactions and/or a staggered component of the g (gyromagnetic) tensor. Such 
terms lead to a gapped quasiparticle spectrum and a nonzero condensate density for all tem- 
peratures so that no phase transition occurs. We contrast this to the effect of crystal field 
anisotropics which are also induced by spin-orbit coupling. In the second part we study the 
field-induced magnetic ordering in TlCuCls on a quantitative level. We show that the usual 
BEC picture does not allow for a good description of the experimental magnetisation data and 
argue that antisymmetric spin interactions and/or a staggered g tensor component are still 
crucial, although both are expected to be tiny in this compound due to crystal symmetries. 
Including this type of interaction we obtain excellent agreement with experimental data. 



Introduction. - It has been argued that the phase transition at a critical applied magnetic 
field in certain spin systems with an excitation gap A as for example integer-spin antiferromag- 
netic chains [1], weakly coupled two-leg ladders [2] or three-dimensional dimer systems [3,4] 
may be regarded as a Bose condensation. In these systems the lowest excited state is a triplet 
of massive bosons. A magnetic field H causes a Zeeman splitting of the triplet with the lowest 
mode crossing the ground state at a critical field He — ^lg\iB- The ground state for H > He 
then becomes a BEC of this low-energy boson. The density n of the boson mode is directly 
related to the magnetisation per site m = gfJ-sn. In principle, spin-gap compounds offer 
therefore the exciting possibility to study BEC in a system where the density n is tunable by 
the external magnetic field which acts as a chemical potential for the bosons. However, as 
the triplet excitations (magnons) interact via a hard-core repulsion a description as a weakly 
interacting Bose gas is only meaningful if the average distance between the magnons / ~ n~^^^ 
is much larger than the s-wave scattering length a which is the characteristic length scale rep- 
resenting the influence of the repulsive potential. This implies that a/l ^ n^^^a <C 1 so that 
the magnons have to be dilute. In this case the well-established gas approximation [5] which 
involves a systematic expansion in terms of the small parameter n^/^a is applicable and even 
the finite temperature properties of the interacting Bose gas can be studied analytically [6] . 
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A dilute magnon gas is realized in the spin dimer system TlCuCla in magnetic fields 
H ^ 6 — 7 T and considerable interest has focused onto this compound within recent years 
[3,4,7-10] . TlCuCls has an excitation gap A w 0.7 meV in zero magnetic field and a bandwidth 
W ^ 6.3 meV [8]. The dimers in this compound are formed by the S = 1/2 spins of the Cu^+ 
ions and superexchange interactions are mediated by the Cl~ ions. The crystal structure can 
be considered as two-leg ladders formed by these dimers. However, inelastic neutron scattering 
(INS) [8] has revealed that the magnons show a considerable dispersion in all spatial directions 
indicating that the interladder interactions are strong. TlCuCls therefore has to be considered 
as a three-dimensional (3D) interacting dimer system. For fields H > He ~ 5.6 T a long- 
range magnetic ordering below some critical temperature has been detected. It has been 
shown that this transition and the overall shape of the magnetisation curves as a function 
of temperature can be qualitatively described as the BEC of magnons [3]. Recently, the 
excitations for H > He have been measured by INS and the lowest mode has been interpreted 
as the gapless Goldstone mode characteristic for a Bose condensed phase [11]. 

In the first part of this letter we will show how to include antisymmetric spin interactions 
or a staggered component of the g tensor into a Hartree-Fock-Popov (HFP) treatment [3,6] of 
an interacting dilute Bose gas. Such terms will in general break the axial U{1) symmetry of 
the system so that a Goldstone mode no longer exists. Furthermore the condensate density 
becomes nonzero for all temperatures so that no phase transition occurs. These findings are 
in agreement with [12] where the same type of anisotropy in a Haldane spin chain has been 
considered. We will discuss the differences to the effect of crystal field anisotropics which 
can also break U{1) symmetry and affect Bose condensation [1]. Our analysis will allow us 
to investigate the effects of such anisotropics on the magnon density in 3D dimer systems 
at finite temperatures quantitatively and we present a detailed study of the field-induced 
magnetic ordering in TlCuCls along these lines in the second part of this letter. 

General scenario. - We want to restrict ourselves here to 3D spin-1/2 systems where 
the spin gap is due to some kind of explicit dimerisation. A useful approach to describe the 
excitations in such systems is the bond operator representation for spins introduced in [13]. 
The starting point is the strong coupling ground state \s) where each dimer at site i is in 
singlet configuration \i,s). It is then natural to introduce operators tj^ which create local 



triplet excitations \i,a) = tlJi,s) with = -\ tt), = \ U) and |z,0) = (| ti) + | it 



))/y/2. Due to hopping between the dimers the three triplet components acquire a dispersion 
flka = ^ko — ctgfJ'sH with a band minimum ^q^a = A — agusH at momentum depending 
on microscopic details of the exchange interactions. The triplets are subject to a hard-core 
constraint which can be taken into account by introducing an infinite on-site repulsion 



with a, P = — , 0, Note that the renormalisation of the triplet dispersion Qka due to the two- 
particle scattering vertex v{k, lu) corresponding to the interaction (QJ can be calculated exactly 
in the dilute limit by a summation of ladder diagrams [14]. For magnetic fields H > He and 
temperatures T < A it is sufficient to take only the lowest triplet mode (a — +) into account. 
In this case only particles near the band minimum at Qq are excited so that the generally 
energy and momentum dependent v{k,Lu) can be replaced by a constant vq — v{qQ,0). With 
the definitions = D,ko — A and = tk+ the Hamiltonian for the lowest mode is given by 




U 



oo 



(1) 



i.a.,(3 



(2) 
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where no — gixsiH — He). Here we want to consider a perturbation linear in i, t'^ 

n' = i7{t,„-tl) (3) 

where 7 is a small parameter. This term respects parity and time reversal symmetry and will 
therefore be non-zero in general if not forbidden by additional crystal symmetries. Clearly 
its effect will be non-perturbative at fields H ~ He because it mixes the singlet ground state 
and the triplet soft mode. Physically such a term can originate from Dzyaloshinsky-Moriya 
(DM) interactions ^ D ■ {Si x Sj) between spins at sites i,j where D is the DM vector. For 
illustration we consider the two-leg spin-ladder shown in fig.^ In the limit J± ^ J the triplet 




Fig. 1 - Ladder with rung couplings J± and J <^ J± along the legs. D is staggered along the x-axis. 
The squares denote some environment preventing inversion centres within each dimer. 

dispersion for this model is given by flka — J± + Jcosfc — agfisH and has a minimum at 
qo = TT. Consider now the additional antisymmetric exchange 

Hdm = ^h^y {S\ySl^ ~ S\^Sly) (4) 

i 

where indices I, r label the two spins within the dimer as shown in fig.H Using the identity [13] 
iS'^''' — {±ta ± tjj — leap^ft'^pt-f) /2 with a = x,y,z this is easily transformed into Hdm = 
iD^^-{—iy{ti^ — t\^)/2. Using spiral instead of cartesian indices we find in momentum 
space Ti-DM ~ ^Dx{t.,r — t^) where t denotes now the lowest triplet mode (a — +) as in |(2Jl. 
This is exactly of the proposed form Q with the triplet operators il^'' acting at the band 
minimum. It can be shown by an exact transformation that the considered DM interaction 
produces an effective staggered field if a uniform magnetic field is applied [15]. A staggered 
magnetic field can also originate directly from a staggered g tensor [12]. Both effects are 
induced by spin-orbit coupling and can contribute to (PJ- Note, however, that the interaction 
((2Jl itself is forbidden by symmetry if there is an inversion center in the middle of each dimer 
irrespective of the origin of this term. 

Next we diagonalize ^ with the perturbation lO included, treating the interaction be- 
tween non-condensed magnons in the one- loop (HFP) approximation. We want to emphasize 
again that taking only these diagrams into account is not an uncontrolled approximation but 
instead the first order in a systematic expansion in the gas parameter n. First, we introduce 
new operators by tk = Ck + iSk g^i] where 77 is a real number. The density of condensed 
magnons per dimer rip is then given by uq = rj^. Ignoring a momentum independent term, the 
Hamiltonian splits into 2 parts H = Tiiin+WbiUn with Hun — K'^nvoV+'^oV'^ ^ l^o'n^l)i'^qQ~'^qo) 
where n denotes the density of non-condensed magnons per dimer and 

Wbiiin - J2iAkcick - ^uicicl^ + h.c.)/2} . (5) 

k 

Here Ak = Cfc— Mo+^ii with the normal self-energy Sn = 2von+2vorj'^ whereas the anomalous 
self-energy is given by E12 = Vo'ij'^. By a Bogoliubov transformation we find the quasiparticle 
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spectrum Ek — {A\ 



S^2)^^^- Additionally we have to demand that Tiun vanishes 



^J■Q'n -l + V (Sii - S12) = . 



(6) 



To stress the differences between BEC with and without TC' we first briefly summarise the 
case 7 = which is the usual BEC scenario for a weakly interacting Bose gas in the HFP 
approximation. In this case we have to distinguish further between the case 1] = 0, i.e., no 
magnons are condensed and the case with nonzero condensate density ry 7^ 0. For 77 = 
the subcondition © is identically fulfilled, the anomalous self-energy S12 vanishes and the 
quasiparticle spectrum is identical to the bare triplet spectrum with an effective chemical 
potential fj, = fiQ — 2vQn. The density of magnons in this phase is therefore given by the 
usual Bose distribution n = {1/N) ^^[exp(/3(efc - - 1]"^ BEC occurs when the chemical 
potential fi vanishes so that the density rif, at the critical point is given by Uc = /io/2wo- At 
temperatures below the critical point (T < Tc) J] will be nonzero and the condition © for 
7 = becomes equivalent to the Hug enholtz- Pines theorem [6] /Ltg = Eii — S12 and guarantees 
the existence of a Goldstone mode, i.e., the quasiparticle spectrum Ek = \/e^ + 2efcnoWo is 
gapless. The number of non-condensed magnons h in this phase is given by 



whereas the number of condensed magnons can be determined from cq. Q which has to be 
solved self-consistently together with iQ. 

Now we want to compare this with the case 7 7^ 0. From eq. ^ it follows immediately 
that this changes the situation qualitatively irrespective of the magnitude of 7 because rj 
must be nonzero for all temperatures so that there will be always condensed magnons and no 
phase transition will occur. Furthermore there is no longer a Hugenholtz-Pincs theorem and 
consequently the quasiparticle spectrum 



is gapped. This is expected because the component of the DM interaction perpendicular to the 
applied field H breaks the U{1) symmetry so that there is no longer a Goldstone mode. Using 
the spectrum (|SJ, the non-condensed magnon density n can be again calculated by eq. (jT)) 
with Efc being replaced by Ck + h\/y/no- The condensed density is determined by eq. © so 
that there are again two equations which have to be solved self-consistently. 

Finally we want to discuss the difference between the perturbation |(2Jl originating from DM 
interactions or a staggered g tensor and single-ion anisotropies ~ D[Sf)'^ + E[{Sf)'^ — )^]- 
(Exchange anisotropies have the same effect as single-ion anisotropies) In terms of triplet 
operators a single-ion anisotropy produces a perturbation bilinear in the triplet operators 
Wpert = "/{tkt-k + h.c.) which rcsults from the _E-term provided that the magnetic field is 
along the z-axis. Here 7 is again a small parameter and we assume 7 > without loss of 
generality. Consequently the condition for the vanishing of the linear terms now becomes 
—fJ.ori — 277/ -f- rj (Sii — S12) = and 7 7^ does not imply ry 7^ as before. Therefore a phase 
transition will still occur and no = for T > Tc- The quasiparticle gap will vanish exactly 
at the critical point, however, it will reopen again below Tc- Additionally, the critical point 
itself will be slightly shifted. Although in both cases the quasiparticle spectrum is gapped for 
no 7^ 0, the gap only weakly depends on the condensate density for the single-ion anisotropy 
case contrary to what was found before (see eq. 






(8) 
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Fig. 2 - Experimental magnetisation curves (symbols) for different magnetic fields H = 
5.9, 6.0, 6.1, •• ■ ,7 T taken from ref. [3] compared to the theoretically calculated (solid lines) with 
A — 0.67 meV where (a) vo — 9.8 meV and (b) vo = 25 meV. In (c) the result with included is 
shown where 7 = 10~^ meV, A = 0.72 meV and vq — 27 meV. (d) Quasiparticle spectrum (|HJ (solid 
line) for = 14 T (no « 2.3 ■ 10"^ at T < 1.5 K) with aU other parameters as in (c). INS data for 
T = 50 mK (circles) and T — 1.5 K (squares) taken from ref. [11] are shown for comparison. 

TlCuCls- - For TlCuCla detailed magnetisation measurements have been performed for 
magnetic fields H > He with typical magnon densities n ^ 10^"^ [3] satisfying the condition 
of diluteness. This compound is therefore particularly suited to test the predictions of the 
BEC theory within the HFP approximation even quantitatively. This is the purpose of the 
remainder of this letter. TlCuCla crystallises in the space group P2i/c with two dimers per 
unit cell formed by the 5 = 1/2 spins of the Cu^+ ions [8]. Because the centre of each dimer 
is an inversion centre, the interaction (O is forbidden by symmetry making, apparently, a 
description as a usual BEC valid. To do so it is not sufficient to use a magnon dispersion 
of the form A + /2m [3]. Although a quadratic dispersion is indeed expected close to 
the band minimum its applicability is restricted to small excitation energies and we find 
that this simplification here is only justified for temperatures T < 1 K which is well below 
the experimental temperature range (}). Matsumoto et al. [7] have successfully applied the 
afore mentioned bond-operator technique to describe the real triplet dispersion in TlCuCls. 
However, we have found that the mean-field treatment of the hard-core constraint used 
in [7] is not sufficient for a quantitative comparison with the magnetisation measurements 
in [3] . We therefore have calculated the renormalisation of the dispersion due to the constraint 
in a systematic way by a summation of ladder diagrams [14]. Doing so we find a considerable 
renormalisation of the superexchange parameters compared to ref. [7]. In particular this 
treatment allows us to calculate the scattering amplitude vq directly and we find vq = 1.6 W = 
9.8 meV [16]. Additionally we note that the gap cannot be determined very accurately from 
the INS data so that we have used He = 5.6 T and g = 2.06 [17] which yield A = gusHc = 0.67 
meV. This value was used as a constraint when calculating the renormalised dispersion. Having 
the renormalised spectrum and wq at hand we can calculate the magnon density and therefore 
the magnetisation as a function of temperature within the standard BEC picture. As we are 
interested in temperatures T <C wo we can treat vq as temperature independent. The result 
is shown in fig. Quite obviously the calculated magnetisation curves do not agree with 
the experimental data. We want to point out that although the jump in the magnetisation at 

(^)The experimentally observed deviation from the universal BEC power law Uc ~ T" with a = 3/2 is 
therefore not astonishing because this exponent is a direct consequence of a quadratic dispersion. 
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the transition point is known to be an artefact of the HFP treatment, this approximation is 
rehable for a dilute Bose gas apart from the small temperature interval \T — Tc\ < n^^^aTc [6] 
where we find n^^^a ^0.1. Especially the huge overestimation of the magnetisation at T = 
is definitely not due to an invalid approximation. 

Taking vo as a fitting parameter we can obtain good agreement with experiment up to 
the minima of the magnetisation curves with vq = 25 meV yielding the best result shown in 
fig. 1213. A rcnormalisation of vq up to this value is by no means implausible as our calculation 
of Vq takes only magnon-magnon interactions into account. Raman spectroscopy, however, 
has revealed that low-energy optical phonons exist in TlCuCls which interact with the spin 
system [9] . As a consequence a reduction of the bare magnon bandwidth Wq due to polaronic 
effects is expected so that t^o, which is basically given by Wq, might be in fact much larger. A 
reduction of the bandwidth by a factor 2-3 is not exceptional and known in polaron physics for 
a long time [18]. However, even when using vq = 1.6 Wq = 25 meV the theoretically calculated 
magnetisations are still about 50% too large at T = 0. More important, the standard BEC 
scenario always predicts a sharp increase of the magnetisation below Tc which qualitatively 
disagrees with experiment where only a slight and smooth upturn is visible. 

In recent ESR measurements [17] a direct singlet-triplet transition has been observed which 
would be usually forbidden by spin conservation. This yields some indication that a small 
perturbation of the form (PJ is still present. The reason might be small static distortions 
which violate the exact inversion symmetry within a dimer so that 7i', although certainly 
tiny, becomes nonzero. We therefore tried to fit the measured magnetisation curves using the 
outlined BEC theory with Ti.' included. Because any magnetic field will then cause a finite 
magnetisation. He is no longer well defined. However, fixing the gap in a range consistent 
with the INS data we can still use He = A/gfiB as a formal definition. With A = 0.72 meV, 
Vq = 27 meV and a DM interaction 7 — 10~^ meV we can obtain excellent agreement with 
experiment as shown in fig.|2t. Most important, even this extremely small perturbation yields 
smooth minima and a slow increase of the magnetisation at temperatures below the minima 
consistent with experiment. If 7 is indeed nonzero we expect a gapped quasiparticle spectrum. 
In fig. OJi we therefore compare the spectrum ((SJ with the INS data taken from ref. [11]. At 
low excitation energies our result is in perfect agreement with these data. The deviations at 
higher energies are expected because our approximation v{k,uj) — v{qQ,0) in eq. (O is then 
no longer justified. The predicted gap A ~ 0.1 meV is only a factor 2 smaller than the lowest 
measured excitation energies making it perhaps accessible in future studies. We also mention 
that in systems where DM interactions are allowed by symmetry they are typically of the order 
7 ~ J {9 ~ where J is the isotropic superexchange [19]. For a DM interaction within a 
TlCuCls dimer this estimate yields 7 ^ 0.2 meV which is two orders of magnitude larger than 
the term considered here supporting our statement that tiny violations of inversion symmetry 
are sufficient. However, a nonzero anisotropy term raises questions about the orientation of 
the DM vector D. li D would be oriented along a specific axis throughout the crystal as in 
our example in fig. ^ it would be possible to restore U{1) symmetry by applying the magnetic 
field along the same axis. In this configuration a sharp phase transition would be visible. 
However, such dependencies on field direction have not been reported for TlCuCls. Another 
possibility would be that small static distortions lead to domains with different orientations 
of £>. In this case a component perpendicular to H could exist for each field direction. 
Provided that indeed no anisotropy with respect to the field direction exists we believe this 
is the most probable scenario taking the tininess of the needed DM term into account. In 
principle it is also allowed by symmetry to construct a DM term from the magnetic field itself, 
i.e. it is possible that the magnetic field induces the anisotropy and therefore determines the 
orientation of D. Finally we want to discuss a single- ion or exchange anisotropy term as an 
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alternative possible perturbation. As already mentioned before a sharp phase transition will 
still occur with such a term included. Additionally, the quasiparticle gap in the condensed 
phase will be almost independent of no and small for a small perturbation. We therefore do 
not expect any qualitative change for the magnetisation curves. Explicit calculations indeed 
show that such a term basically yields a small shift of the critical temperature leaving the 
shape of the magnetisation curves otherwise unchanged (data not shown). 

Conclusions. - In summary, we have studied how a term H' caused by antisymmetric 
spin interactions and/or a staggered g tensor affects BEC in spin-gap systems. This interac- 
tion directly mixes the singlet with the triplet soft mode so that its effect is non-perturbative 
at magnetic fields iJ ~ iJc- Treating the magnon-magnon interaction in the one- loop approx- 
imation we have found a gapped quasiparticle spectrum compared to the gapless spectrum 
without Ti.' and a nonzero condensate density for all temperatures consistent with findings in 
ref. [12] for Haldane spin chains. Differences to the case of crystal field anisotropics have also 
been discussed. We have then performed a detailed quantitative study of the field induced 
magnetisation process in TlCuCls and found that the scattering amplitude vo seems to be 
much larger than expected from magnon-magnon interactions alone. We see this as a fur- 
ther confirmation that spin-phonon coupling is important in this compound. Finally we have 
pointed out that even in TlCuCls where H' is expected to be tiny due to symmetry it seems 
to remain crucial to obtain a quantitative correct description of the experimental data. 
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